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In a quantum cosmological model consisting of a Euclidean region and a Lorentzian region, Hartle- 
Hawking's no-bounary wave function, and Linde's wave function and Vilenkin's tunneling wave 
function are briefly described and compared with each other. We put a particular emphasis on 
semiclassical gravity from quantum cosmology and compare it with the conventional quantum field 
theory in curved spacetimes. Finally, we discuss the recent debate on catastrophic particle produc- 
tion in the tunneling universe between Rubakov and Vilenkin within the semiclassical gravity. 
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I. INTRODUCTION 

It is widely known that quantum gravity effects are 
important at the Planck scale of length, time and mass 
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It is also well accepted that at such Planck scales quan- 
tum gravity theory should replace the role of classical 
gravity theory. Though at present any viable theory of 
quantum gravity is not known free from all conceptual 
and technical problems, quantum gravity may be denoted 
symbolically as 



(2) 



It is within quantum gravity with both gravity and mat- 
ter quantized that the unitarity problem and renormal- 
ization problem, quantum black holes, and quantum cos- 
mology should be properly addressed. 

In a regime where gravity becomes classical but mat- 
ter still maintains quantum nature, one has semiclassical 
gravity theory 
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Gravity may undergo, for instance, a decoherence mech- 
anism to make a quantum-to-classical transition. Quan- 
tum field theory in curved spacetimes or the semiclassical 
theory derived from quantum gravity may be regarded 
as semiclassical gravity theory. This semiclassical grav- 
ity theory has been widely used to investigate important 
phenomena such as particle and/or entropy production. 



black hole information, and some of cosmological prob- 
lems. Finally, when both gravity and matters become 
classical, they are governed by classical gravity theory 
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Even with a consistent theory of quantum gravity 
known, one should also understand the mechanism lead- 
ing quantum gravity 10) to semiclassical gravity Q and 
finally to classical gravity l@J. The main aim of this talk 
is to follow this scheme, discuss and clarify some issues 
in this direction P, As we do not know the con- 
sistent quantum gravity, we need to adopt some quan- 
tum gravity theory that incorporates quantum gravity 
effects. The Wheeler-DeWitt (WDW) equation, for in- 
stance, may be such a quantum gravity theory where 
both gravity variables from metric tensors and matters 
are quantized. Quantum cosmology, where cosmological 
models are quantized a la the WDW equation, can be a 
good arena to test the above scheme of a transition of 
quantum gravity down to classical gravity. 

Quantum cosmology consists of a quantum dynamical 
law, boundary (initial) conditions and prediction (inter- 
pretation) of quantum states Q . The WDW equation as 
a quantum dynamical law determines the wave functions 
of the universe. There should be some prescription to 
select a unique wave function of the universe. A bound- 
ary condition selects the unique or appropriate state of 
the universe. Then the selected wave function carries 
all information of the history of the universe. Three 
leading proposals to prescribe boundary conditions were 
advanced: Hartle-Hawking's no-boundary wave function 
0, Linde's wave function |^, and Vilenkin's tunneling 
wave function |^. We derive the semiclassical gravity 
from the WDW equation and discuss particle creation in 
a tunneling universe. Recently there has been a debate 
on catastrophic particle creation from a tunneling uni- 
verse between Rubakov's group and Vilenkin's group 
0. We show that an adiabatic vacuum state minimizes 
particle creation but its squeezed states necessarily lead 
to catastrophic particle creation. 
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II. WAVE FUNCTIONS OF THE UNIVERSE 

As a simple model for quantum cosmology, we consider 
the Friedmann-Robertson- Walker (FRW) universe with 
the action 



167r 



R-2K 



+ ^ I d^xVhK, (5) 
57r 



where A is a cosmological constant. The surface term 
is introduced to yield the correct equation of motion for 
a closed universe. In the Arnowitt-Deser-Misner (ADM) 
formulation with the metric 



ds^ = -N^{t)df + a^{t)d^l, 
the action leads to the Hamiltonian constraint 
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pf A 



= 0, 



where the conjugate momentum is given by 

S-TTTOpa da 
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The quantum theory of the universe is given by the 
Hamiltonian constraint a la the Dirac quantization, 
which is nothing but the WDW equation 



2m^ da? 



+ Vg{a) ^{a)^Q, 



with an effective gravitational potential 
' 8 V 3 
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The WDW equation for pure gravity of the FRW uni- 
verse now has the form of one-dimensional Schrodinger 
equation with the potential H10(l and zero energy. No ex- 
act solutions are known for this apparently simple equa- 
tion. More significantly, the boundary conditions for 
quantum cosmology may differ from those for quantum 
mechanics. However, the characteristic behaviors of the 
wave functions of leading proposals can be well exhibited 
by the asymptotic solutions For that purpose, we 
introduce the Liouville- Green transformation |lO| 
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and write the gravitational wave function in the form 

-d^ 
-drj'^ 



+ T] + A{r]) ^(77) = 
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where 
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Here at is the turning point, a < at = (A/S)""^/^, from the 
side of a de Sitter region. The Liouville- Green transfor- 
mation has an advantage that it can readily be extended 
to classically forbidden regions without singular behavior 
at turning points. The integral (|11|1 can be done exactly 
to yield 
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We can find asymptotic solutions to Eq. ^ at least 
in two asymptotic regions, a de Sitter region and a clas- 
sically forbidden region. In the first case of the de Sitter 
region, (77 > 1,A = ~/16r]^ 0), Eq. (O has the 
asymptotic solutions 
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where Ai and Bi are Airy functions and (5 is a small 
correction due to A. In the second case of the classically 
forbidden region of tunneling universe, we analytically 
continue Eq. (|ll|l to introduce another variable 
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and write the gravitational field equation as 

d^ 
VdC 



C + A(C) *(C) = o, 



where 

3(dW) 
^ ' 4 ^dC/daY 

Now Eq. Iini) yields 
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The asymptotic solutions to Eq. (|17|l in the tunneling 
region are 



A^{c-m) 
B^{c-m) 
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Using the asymptotic solutions (|15|l and (|2()|l . we 
can prescribe the three leading proposals for the wave 
functions of the universe. First, Hartle-Hawking's no- 
boundary wave function in the classically forbidden re- 
gion is given by 



^HH{a)^CAi{C-5{C)), 



(21) 
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and has the asymptotic form 



1 \l/4 1 



(22) 



Here Sg is the gravitational instanton 



Sg{a)^£{2mlVg(a)f'^fe/'. (23) 

In the de Sitter region we analyticaUy continue C to 77 to 
obtain 



^'ffH(a) = CAi{-'q- S{ri)). 



(24) 



In fact, the wave function H24|l is an analytical continu- 
ation of H21|l through the turning point without having 
any singularity. The Hartle-Hawking wave function has 
the asymptotic form 



*HH(a) = 
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sin(-5g(a) + -), (25) 



where 



Sg{a) = / {~2mlVg{a)) ' - —77' 



Second, Linde's wave function is prescribed by 

and has the asymptotic form 
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where Sg is the gravitational instanton H23|l . The wave 
function can be continued to the de Sitter region 



with the asymptotic form 

1 \l/4 



4'L(a) 



(29) 



(30) 



Finally, Vilenkin's tunneling wave function is given by 

^v{a) = C[Ai{C - ^(0) - iBliC - 6{0)], (31) 
with the asymptotic form 

vl/^(a) = c( [ie-*^»(") + . (32) 

V 1 — jQ-^ ' L2 J 

In the de Sitter region it becomes a purely outgoing wave 
function 



(33) 



A few comments are in order. First, in the tun- 
neling region, Hartle-Hawking's, Linde's, and Vilenkin's 
wave functions have the asymptotic forms, "^hh ~ 
^-s,{a)lh^ ^ ^+S,/h^ g^^^ ^ ^+S,/h^^-S,/h/2, re- 
spectively. As 5*3(0) is a decreasing function of a, Hartle- 
Hawking's wave function increases from a minimum value 
whereas Linde's and Vilenkin's wave functions decrease 
from a maximum values at a = 0. When 5*^ is large, 
Linde's wave function and Vilenkin's wave function are 
dominated approximately by the same asymptotic form 
gSg/h^ These forms are important in deriving semiclassi- 
cal gravity in the next sections. Second, in the de Sitter 
region Hartle-Hawking's wave function and Linde's wave 
functions have equal amplitudes (probabilities) for out- 
going (expanding) and incoming (collapsing) components 
whereas Vilenkin's wave function has a purely outgoing 
component. 



III. SEMICLASSICAL GRAVITY 

The idea of deriving semiclassical gravity from quan- 
tum gravity is rooted on the correspondence in certain 
circumstances between quantum theory and classical the- 
ory. In fact, the Schrodinger equation 
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V{x) ijjix) = EiPix), 



with a WKB solution of the form 



may lead, in the asymptotic limit (l/h 
Hamilton- Jacobi equation 



1 
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iVS'f + V{x) ^E. 



With the identification p — VS" and B 
becomes nothing but the Hamiltonian 

B = + V{x). 

Zm 



(34) 
(35) 
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The same idea may be applied to quantum cosmology. 
The WDW equation for gravity ha coupled with a matter 
field (j) takes in general the form 



+ mlV{ha) + Brnin^, ^, ha) *(/la, ^) ^ 



2m^ 



Here ml = l/G, and 



(38) 

— ^1 ^, ^..^ , - GabS^ /ShaShb is the Laplace- 
Bcltrami operator on the superspace ha of three geome- 
tries modulo diffeomorphisms, and V{ha) is a potential 
from the three-curvature. Contrary to ordinary quantum 
theory, there are two asymptotic parameters, mp — > oo 
and 1/h oo. A simple analogy is to write the wave 
functional as j3| 



vI/(/i„(/.)=e^^°(''")x(0,/ia). 



(39) 
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In the limit of — > oo, one obtains the Einstein- 
Hamihon-Jacobi equation 



1 



2m2 



iS/Sof + mlV{ha) = 0, 



and the time-dependent Schrodinger equation 

i^^X = -ffm(7r0,0, ha)x- 



(40) 



(41) 



Here a cosmological time is introduced along each clas- 
sical trajectory determined by the Einstein-Hamilton- 
Jacobi equation as 
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However, this simple idea has a limitation that Eq. H40|l 
is a pure gravity (vacuum solution) without the backre- 
action of matter. 

To overcome this shortcoming, we may employ the 
Born-Oppcnheimer idea of separating gravity from mat- 
ter, and the de Broglie-Bohm interpretation of emergence 
of classical system with a quantum potential 0. That 
is, the wave function is written as 



^'(/la,0) = 1p{ha)^i(f>,ha), 



(43) 



where ip{ha) is a function of the slow variable ha with 
mass scale nip and ^{(p, ha) is a function of the fast vari- 
able (f> with mass scale m. Then one may obtain from 
the WDW equation the gravity part equation 



D' + mlV+{H„ 



and the matter part equation 

h^ 1 



2ml 



(52) V = 0, (44) 



(£)2 _ (£)2))$ ^ 0. 



(45) 
(46) 



Here D and D are covariant derivatives 



and dropping tildes, and adopting the de Broglie-Bohm 
interpretation 

,Piha) ^ Fiha)ei''^''-\ (51) 
one finds the Einstein-Hamilton-Jacobi equation 
1 



2m2 



(52) 



and the quantum field equation for the matter field 



ih-^ = [H„, - (i?„)]$ = 0. 



(53) 



The cosmological time is defined with respect to the grav- 
itational action with the backreaction from the matter 
field 



1. J-VS-V. 

dt mi 



(54) 



In conclusion, the semiclassical gravity from quantum 
gravity takes the form 



mi 



^h-m ^ 



Hm + H2 



|1>) =0, 



(55) 



where and are quantum corrections. For in- 
stance, the semiclassical gravity for the FRW cosmology 
is given by 
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(56) 



where 



Un 



1 d{ad)/dT 



(57) 



Da — ^ i-^a. Da — "r; iAa, 

Oha Oha 



with an effective gauge potential 



Aa{ha) = -i 



and {H„i) and (I?^) are backreactions 



($1$) ' ^ ' ($1$) 

After removing the geometric phases 
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(48) 
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IV. SEMICLASSICAL GRAVITY FOR 
TUNNELING UNIVERSE 

From now on we shall confine our attention to the FRW 
universe with a minimal scalar field. The action for a 
minimal massive scalar field 



d^^{-9)'^^l [-g^-'d^^d^cf) - m202 



(58) 



after the Fourier-decomposition, leads to the canonical 
Hamiltonian of the form 
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Here a denotes the collective notation for the Fourier or 
spherical harmonic modes 



(for a spatially flat FRW), 
2), (for a closed FRW), 



(60) 



and TTa = a^4>a is the canonical momentum. Thus the 
Hamiltonian is equivalent to an infinite sum of time- 
dependent oscillators. The Hamiltonian constraint for 
gravity plus matter system now becomes 



Hm — 0, 



(61) 



and leads to the WDW equation 



3a 



(62) 

Without complication of a backreaction of (/> and loss 
of generality, we assume that only a few modes be present 
in Eq. H62|l . We then define the cosmological time 142(1 
along the classical trajectory of an expanding universe 
with the gravity part of wave function 



V'cxp(a) =F(a)e*^«(") 



as 



d_ 

dt 



dSg d 



STrrripa da da 



Then the field equation becomes 



d 



ih—^{(j)a,a{t)) ^ Hm{(l)a,a{t))^{(j)a,a{t)). 



(63) 
(64) 

(65) 



This expanding wave function can be analytically con- 
tinued along a semicircle in the upper half of a complex 
plane to the gravity wave function in a tunneling region 



V'. 



HH 



F{a)e 



(66) 



Then the cosmological time in the tunneling universe be- 
comes 



d_ 



dS„ 



Sirmpa da ' 



(67) 



implying an imaginary time t = it. The field equation 
now takes the form 

n^$(0a,a(T)) = i7„,(0„,a(T))$((/)„,a(T)). (68) 

OT 

Note that the quantization rule in the tunneling universe 
(Euclidean geometry) changes to 



(69) 



In the coordinate representation, tTq, — hd/d(f)a, the field 
equation in the tunneling universe becomes 



d 

h—^4>a,a{T)) 

OT 
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2a3, 



2 i 2 



$(0„,a(T)). 

(70) 



In passing, we note that the time-reversal operation of 
field equations can be naturally understood in quantum 
cosmology. The wave function for the collapsing universe 



(71) 



with respect to the cosmological time leads to an- 

other field equation 

-ih—^{(t)a.,a{t)) = Hm{K,a{t))^(t)c.,a{t)). (72) 

In fact, Eq. |(721) is the time-reversal, $(i) = ^*{-t), of 
the field equation l(65|l from the expanding universe. Sim- 
ilarly, in the tunneling universe the gravity wave function 



V' = F(a)e+*^''("), (73) 
leads to another field equation 
d ~ 

-ft $(^^,a(r)) = H„^{(t>o.,a{T))<^>{4)o,,a{T)). (74) 

OT 

Note that the field Equation (|74|l . being real, is also the 
time-reversal of Eq. (|68|l . However, this equation is nec- 
essary to equip the quantum field theory with an inner 
product 



(($2 (</'o,t)|0|$i ((/.„, r)) = j d0,$2(</'a,T)O$i((/.„,r). 

(75) 

It should be remarked that, with the time convention 
(|67|l . Hartle-Hawking's wave function leads to Eq. H68|l 
whereas Linde's wave function leads to Eq. (|71l) . These 
two equations complement each other with respect to the 
inner product (|75|l . 



V. QUANTUM FIELD THEORY: LORENTZIAN 
VS. TUNNELING UNIVERSES 

In this second part we study and compare the semi- 
classical gravity in a tunneling universe and a pure 
Lorentzian universe. We first introduce a scheme in the 
Lorentzian universe to find the Fock space of the min- 
imal massive scalar field in Sec. HI, and then extend 
the scheme to the tunneling universe to obtain the cor- 
responding Fock space. 



A. Fock Space in Lorentzian Geometry 

The Hamiltonian H59I) of the minimal scalar field con- 
sists of time-dependent oscillators that depend on time 
through a. The invariant method |llj | turns out a pow- 
erful tool in finding various quantum states for such a 
time-dependent quadratic Hamiltonian. Following Ref. 
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|l2j |. we introduce the time-dependent annihilation oper- 
ator 



[ip*Jt)iT^-aHt)rjt)<t>^], (76) 



and the creation operator the Hermitian conju- 

gate of Aa(t), and require them to satisfy the quantum 
Liouville-von Neumann equation 



ih—Ao,{t) + = 0. 



(77) 



Here an overdot denoting a derivative with t. Then the 
auxihary field satisfies the classical equation for each 
mode, 



LPa +3-(Pq +UJ„Lpa = 0. 

a 



(78) 



In fact, these operators satisfy the equal-time commuta- 
tor 

K(t),4(0] =<5./3, (79) 
when Lpa is restricted by the Wronskian condition 



(80) 



The number state of the number operator, Na{t) = 



(81) 



is an exact solution to Eq. H65(l . Once we select a solu- 
tion ifao satisfying both Eqs. ifT^ and l|Hn)l. we find the 
gaussian and its excited states as 



2"°(n„)!V2^p„o 



1/2 



-i0Qo(nc,o + l/2) 



exp 



(82) 



where (pao = Paoe~^^°'° . The most general gaussian state 
and its excited states can be obtained from a more gen- 
eral complex solution of the form 



(83) 



where and Va are time-independent parameters de- 
fined as 



jia — coshrQ, 



sinh Ta ■ 



(84) 



In fact, the time-dependent annihilation and creation op- 
erators in Eq. I|76|l defined in terms of Lpao and ipav are 
related through the Bogoliubov transformation 

Aau{t) = (cosh r„ )i„o (i) - (e**° sinh r„)ilo(t), 

iL(0 = (cosh r„ )ij;,o (t) - (e^*''° sinh r„)i„o(t). 

(85) 

The zz-dependent number state that is obtained by re- 
placing Lpan with Lpai, in Eq. H82|) is a squeezed state of 
the number state (13)] . 



B. Fock Space in Euclidean Geometry 

In the Euclidean geometry of the tunneling universe, 
the non-unitary evolution equation (|68ll has also a pair 
of the invariant operators of the form (14| 

ia+(T) = ---^[(^„+(T)7r„ - a3(T)93„+(T)(^„], 
ia_(r) = -^=\lfia-{T)f^a- 0?{'r)'fa-{'r)(i)a], (86) 

where an overdot denotes now a derivative with respect to 
r and ^Pa± are two independent solutions of the classical 
field equation 



a 



(87) 



These operators satisfy the quantum Liouville-von Neu- 
mann equation in the Euclidean geometry 

?i^ia±(r) + [ia±(r),7?„(r)] = 0. 
With the Wronskian condition 

imposed, the equal time commutation relations 

[A^-{t),Ap+{t)]=5o,P 



(89) 



(90) 



are satisfied. Note then that Aa- (r) is an analytical con- 
tinuation of the annihilation operator A^ (t) and Aq,+ (t) 
is another continuation of the creation operator Al^ {t) of 
the Lorentzian geometry. 

The zero eigenvalue ket-state of Aa-{T), 



A„_(T)|0„,r) =0, 
leads to a gaussian wave function 
1 \i/2 



$ao(0a,T) 



exp 



(91) 



(92) 



As Aa±{T) map one exact solution into another, we ob- 
tain an excited ket-state given by 



1 



na,T 



whose wave function is [If 



:(i.+ (r))"°|Oa,T), 



(93) 



1 



1/2 



2"° {nay.^y2TTh(pa+(pa-' 
(n„-|-l/2)/2_ ^ ^ ra3^^ 



H„ 



exp 



(94) 



One can also show through direct calculation that the 
wave functions H94|) indeed satisfy Eq. (|68|) . 
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The bra-state, which is needed for the inner product 
of quantum states in Eq. H75|l . may be defined by the 
equation 



The vacuum state at r can be expressed as a superposi- 
tion of the number states at tq. The probability is now 
given by projecting into number states at tq as 



and 



((0„,r|ii+(r)=0, 
((n„,rH^(0„,r|(iL(r))' 



(95) 



(96) 



((Oa, r|nQ, ro)((nQ, tqIOq, t) 



(103) 



Here the Hermitian conjugates of ^]j_|_(t) are obtained 

by substituting —ifa for tTq in Eq. In fact, AI±{t) 

satisfy the Liouville-von Neumann equation for the field 
equation 

-ft|:il±(r) + [Al^{r),H^{T)] = 0. (97) 

The wave function for the number bra-state now takes 
the form lla 



It should be remarked that ((Oq, T|nQ, tq) is not a mere 
complex conjugate of {{ua, to|Oq, t), which is true in the 
Lorentzian case. The nonvanishing values are given by 

m.rM = -i,(e2^)"^(-|;-)'nio5) 

Using , we find the mean number of created particles 



(2n, 



ri„(0a,'r) = ( 



1 



2"° {na)ly/2TThipa+ipa^ 
(n„ + l/2)/2 



1/2 



Hn 





r aVa+ 2] 


) exp 





The inner product satisfies the orthonormality 



((ma,T\na,T) = / df?!)^ ^.^^ ((/)„, t) (0a , t) =(5m„n„- 

(99) 



(106) 

These two results are consistent with each other. The 
amplitude squared of Eq. (|104|1 gives a result different 
from H102|l and (|106|l . Thus a complete quantum theory 
in the Euclidean geometry requires both the ket- and bra- 
states, which are not complex conjugates of each other. 
(98) The wave function should be continuous across the 
boundary between the Lorentzian and Euclidean geome- 
tries, $a(to) = *i'Q(''o)- The continuity of the auxiliary 
fields, ifaito) = ifaiTo), implies the continuity of wave 
function. First, we select the adiabatic solution to Eq. 
(|S7|l of the form 



C. Mean Number of Created Particles in 
Tunneling Universe 

As Aa±{T) depend on r through a(T), there is a dy- 
namical squeezing of an initial state. It is the pair Aa- (r) 
and AJj_|_(t) that satisfy the Bogoliubov transformation 
between operators at two different times: 

Aa-ir) = flair, Ta)Aa-{Ta) + I^a{T,To)Al^^{To), 
^i+i'^) = l^a{T,TQ)Al^{To) + i>a{T,TQ)Aa-{TQ). 

(100) 

Though [Aq_(t), AJj_|_(t)] 7^ 1, we can show by direct 
calculation that the Bogoliubov relation holds 



(107) 



where 



-Da±(ro) = \/2a^{To)LOaiTo)ipa±iTo)- (108) 



Then the wave function for the adiabatic vacuum state 
is 



\ 1/4 




j exp 


[ 2h M 



(109) 



Note that Eq. ifTHilll can be obtained through an ana- 
lytical continuation of the adiabatic vacuum state in the 
Lorentzian geometry |l7j | . Now, the mean number of cre- 
ated particles can be calculated by expanding the wave 
function with respect to the number states at the moment 
of entering the tunneling universe 



j-^a^a ^a^a 1- 



(101) 



The initial vacuum state at tq contains particles at r of 
the amount 

{{Oa, To\Al^^{T)Aa-{T)\Oa, To) = i>Q (t, Tq ) J^q (t, Tq ) . 



*Qo(0a,T) = ^C„^$„^(0Q,ro) = y^^Cn^^n^{(j>a,to). 

We thus see that the adiabatic vacuum state minimizes 
particle creation agreeing with the result by Vilenkin et 
(102) al 0. 
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Second, there is one complex-parameter squeezed gaus- 
sian state which can be obtained from a superposition of 
adiabatic solutions 



(cosh TcC^ ^-0 Da± {tq ) 
+6^^^° sinhr„e^'''->°i?„^(To)). (Ill) 



This leads to the squeezed gaussian wave function 



^•ao(</'Q,T) ^TVa+exp 



2h 



(112) 



where Ma+ is a normalization factor and 

c-iSo smhraDa- ~ coshro-e '^^^o Da+ 



sinh rQ,Dct_ + coshr^e 



D 



a+ 



Note that 
Re(B„,) > 



tanhrQ, > e 



-2 

■I T| 



(113) 



(114) 



Thus for a sufhciently long duration of tunneling, even a 
small squeezing parameter eventually leads to an ex- 
ponentially growing state at the end of tunneling. There- 
fore, the squeezed state satisfying the inequality 1)114(1 
forms a one-parameter family of gaussian wave function 
that leads to catastrophic particle production. 



ical constant, we have prescribed three leading proposals 
of Hartle-Hawking's no-boundary wave fucntion, Linde's 
wave function, and Vilenkin's tunneling wave function. 
We also derived in a general context the semiclassical 
gravity from quantum gravity based on the Wheeler- 
DeWitt equation. Further, the difference of semiclassical 
gravity was discussed for the three leading wave func- 
tions. 

Particle creation has been studied for a minimal mas- 
sive scalar field in a tunneling universe from the semi- 
classical gravity point of view. Using the invariant 
method we found the Fock space of the scalar field in 
the Lorentzian and Euclidean geometries of the tunnel- 
ing universe. It was found that the adiabatic vacuum 
state minimizes particle creation in the tunneling uni- 
verse whereas its one-parameter squeezed gaussian states 
lead to catastrophic particle creation for a sufficiently 
long duration of tunneling regime. Thus the debate of 
catastrophic particle creation reduces to the selection of 
gaussian states. 

It would be interesting to compare the result of this pa- 
per with another cosmological model where a tunneling 
region is sandwiched between two asymptotically static 
Lorentzian regions. In this new model the initial vac- 
uum state evolves through the tunneling region into a 
squeezed gaussian state always with catastrophic parti- 
cle creation whereas the evolution without the tunneling 
region suppresses particle creation [l5|. 
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